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We show that for small values of the chemical potential the curvature of the phase transition 
line can be deduced from an analysis of scaling properties of the chiral condensate and its 
susceptibilities. We make use of a recent analysis of the magnetic equation of state in (2+1)- 
flavor QCD where a connection between the QCD parameters and the universal scaling fields 
could be established. The remaining dependence of the reduced temperature on the chemi- 
cal potential can be fixed by an analysis of a mixed susceptibility, obtained from a derivative 
with respect to quark mass and chemical potential. We extract this dependence which de- 
scribes the curvature of the phase transition line, at two values of the cut-off, aT — 1/4 and 
1 /8. We find that cut-off effects are small for the curvature parameter and determine the tran- 
sition line in the chiral limit to leading order in the light quark chemical potential. We obtain 
TM/T,{0) =. 1 -0.00656(66) (mb/T^)^ + ^(Mb)- 
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Figure 1: The curvature of the transition Une in the (r,jUB)-diagram for infinitely heavy quarks and - as the 
result of this paper - in the chiral limit. For comparison we also show the freeze-out data from heavy ion 
experiments, together with a parameterization for the freeze-out line from [ pi] ] . Units are normalized to the 
transition temperature at = (To). 



1. Introduction and Summary 

Extending lattice QCD calculations to non-zero baryon-chemical potential or, equivalently, 
to non-zero net baryon number density is known to be difficult in general. However, important 
information on the QCD phase diagram can be deduced for small values of the chemical potential 
by using well established numerical techniques such as reweighting analytic continuation ^] 
or Taylor expansion [Q, |5|]. At non-zero values of the chemical potential a phase boundary in the 
temperature and chemical potential parameter space of QCD is well defined only in the heavy quark 
limit or for vanishing quark masses. In the former case the phase transition line corresponds to the 
first order deconfinement transition in the pure gauge theory. At infinite values of the quark mass 
this transition is independent of the chemical potential and defines a straight line in the T-pL plane. 
For a large range of quark mass values the transition line is not unique. It characterizes a region 
of (rapid) crossover in thermodynamic quantities and a pseudo-critical temperature extracted from 
these observables may differ somewhat, depending on the observable that is used. In the chiral 
limit, however, the transition line is again well defined. For sufficiently large strange quark mass it 
defines a line of second order phase transitions in the universality class of three dimensional 0(4) 
symmetric spin models [^]. 

In a recent work [Q] we have shown that the curvature of the phase transition line in the chiral 
limit can be obtained from an analysis of the universal scaling properties of a certain mixed suscep- 
tibility which is defined by the leading order Taylor expansion coefficient of the chiral condensate 
with respect to the light quark chemical potential. Numerical calculations have been performed for 
(2-i-l)-fiavor QCD keeping the heavier strange quark mass close to its physical value and decreasing 
the two degenerate light quark masses towards the massless limit. We will make use of a recent 
scaling analysis [^, ^ of the chiral order parameter performed with an improved staggered fermion 
action. This study showed that the chiral order parameter is well described by a universal scaling 
function characteristic for a three dimensional, 0{N) universality class. As a result for the critical 
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line in the chiral limit we find 

T{^b)/T{Q) = l-0.00656(66)(AtB/r)2 + ^((AiB/rf) • (1.1) 

This curvature is about a factor two larger than the reweighting results obtained in (2+l)-flavor 
QCD [[T]]. It is however consistent with results obtained in calculations with imaginary chemical 
potentials. In fact it lies in between the 2-flavor ^ and 3-flavor simulations performed with 
the standard staggered fermion formulation and also is consistent with results reported from (2+1)- 
flavor simulations with imaginary chemical potential performed with the action used also in this 
study (p4-action) [|^]. The result is most relevant for comparison with the experimentally deter- 
mined freeze-out curve as shown in Fig. [l[ The parameterization found in [11| has a curvature 



that is about a factor 3-4 larger which suggests that (for sufficiently large values of the chemical 
potential) the freeze-out line does not follow the critical line. In the following two Sections we 
briefly review the scaling analysis of the order parameter and the mixed susceptibility. 



2. Magnetic equation of state 

In general we can separate two kinds of contributions to the free energy density, a part {fs) that 
will generate singularities in higher order derivatives of the partition function and a regular part 
{fr), we have 

f{T,mum,,lici,lis) = L{T,mi,m„Hg,n,)+fr{T,mi,m„Hy,iJ.,) . (2.1) 

In addition to the temperature T, light (m/) and strange (m^) quark masses we also allow for a 
dependence of the free energy density on the quark chemical potentials. Close to the chiral phase 
transition temperature at vanishing chemical potential the singular part /, will give rise to universal 
scaling properties of response functions. This has been exploited to analyze basic universal features 
of the QCD phase diagram close to criticality [|l2|]. Although fg depends on many parameters of 
the QCD Lagrangian, the universal behavior can be expressed in terms of only two relevant scaling 
variables t and h, that control deviations from criticality at {t,h) = (0,0). To leading order the 
scaling variable h depends only on parameters that break chiral symmetry in the light quark sector, 
while t depends on all other couplings. In particular, t will depend on the light quark chemical 
potential while h remains unaffected by this in leading order, 

h = , (2.2) 

ho Ms 

where Tc is the phase transition temperature in the chiral limit and Iq, ho are non-universal scale 
parameters (as is Tc). While the combination zo = h^J^^ /to is unique for a given theory, the values of 
tQ and /lo will change under rescaling of the order parameter ||8|]. Just like the transition temperature 
Tc also to and Hq are cut-off dependent and will need to be extrapolated to the continuum limit. 

The singular part of the free energy, fs, is a homogeneous function of its arguments. This can 
be used to rewrite it in terms of the scaling variable z = t /h^l^^ as 

f,{t,h) = h'+'I'Uz, 1) ^ h'+"'m . (2.3) 
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where j8, 5 are critical exponents of the three dimensional 0{N) universality class [13]. 

The universal critical behavior of the order parameter, M ~ df/dmi, is controlled by a scaling 
function fc that arises from the singular part of the free energy density after taking a derivative 
with respect to the light quark mass, 

M{t,h) = h'/^fciz) . (2.4) 

The scaling function /g(z) is well-known for the 0(2) and 0(4) universality classes through studies 
of three dimensional spin models [pl|]. This so-called magnetic equation of state, Eq. (2.4), has 



been analyzed recently for (2-i-l)-flavor QCD using an improved staggered fermion formulation 
(p4-action) on lattices with temporal extent A^^ = 4, 8 [0, |8|] and light quark masses as small as 
mi/nis = 1/80, which corresponds to a pion mass that is about half its physical value. It could be 
shown that the chiral order parameter can be mapped onto a universal 0(2) scaling curve ^ and the 
scale parameters to, ho, could be extracted. The scaling analysis has been performed for two 
order parameters, that are multiplicatively renormalized by multiplying the chiral condensate with 
the strange quark mass, but differ in handling additive divergences, linear in the quark mass^, 

A^^^((w)/-^(w).s) • (2.5) 

All resulting fit parameters from fits with and without a regular contribution are summarized in 
Table ^. Note that for = S only fits including a regular contribution have been possible, since 
the smallest available mass in this case has been mi/ifis = I /20. 



3. Curvature of the critical line 

At leading order the light quark chemical potential only enters the reduced temperature t, as 



introduced in Eq. (2.2). Also at non-vanishing values of the quark chemical potential the phase 



transition point is located at f = 0. The variation of the transition temperature with chemical poten- 



tial therefore is parameterized in terms of the constant Kg introduced in Eq. (|2.2D, 



To determine the chiral phase transition line in the T-jj. plane we thus need to determine the propor- 
tionality constant Kg. This is, in fact, the only left over free parameter in universal scaling functions 
that needs to be determined. 

The constant Kq can be determined by analyzing the dependence of the chiral condensate on 
the light quark chemical potential. To extract information about the dependence of the scaling 



' As we are working with staggered fermions 0(2) is expected to be the relevant universality class at finite lattice 
spacing. We find, however, that fits to the 0(4)-symmetric model work equally well. 

^At finite value s of the cut-off these terms are, of course, finite and may be viewed as a specific contribution to the 
regular part that will not alter the scaling properties for sufficiently small values of the quark mass. 
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Mi 


to 


h) 


Tc{Q) [MeV] 


zo 


fit using the scaling term only 


4 


Mh 
M 


0.0037(2) 
0.0048(5) 


0.0022(3) 
0.0048(2) 


194.5(4) 
195.6(4) 


6.8(5) 
8.5(8) 


fit using scaling and regu 


ar terms 


4 

8 


Mh 
M 

Mh 
M 


0.00407(9) 
0.00401(9) 
0.00271(21) 
0.00302(22) 


0.00295(22) 
0.00271(20) 
0.00048(9) 
0.00059(10) 


194.9(2) 
194.8(2) 
174.1(8) 
175.1(8) 


7.5(3) 
7.2(3) 
3.8(5) 
3.8(4) 



Table 1: Scale parameters determined from the scaling fits on lattices of temporal extent Ni ^ A and 8. 
The last column gives zq = h^J^^ /t^. We give the results for parameters entering the definition of scahng 
functions for Mf, and the subtracted order parameter M as defined in Eq. (2.5). Only the former has been 
used in the analysis of the mixed susceptibilities. Note that fits including regular terms, give consistent 
determinations of the parameters of the scaling functions determined from and M, respectively. 



variable t on Kq it suffices to consider the leading order Taylor expansion coefficient of the chiral 
condensate, 



where 



(3.3) 



T di^ijTY dmi/T ■ 

The mixed susceptibility x^^q is proportional to the leading order coefficient of the Taylor expan- 
sion of the chiral condensate, which has been introduced in [|T^, 16]. It may also be viewed as the 



quark mass derivative of the light quark number susceptibility {Xq). 

In the vicinity of the critical point the mixed susceptibility can be expressed in terms of the 
scaling function /^(z) =d/G(z)/dz, 

^ = ^-^h-^'-^)l^'Uz) . (3.4) 

The scaling function /^(z) is easily obtained from fciz) by using the implicit parameterization for 
the latter given in Ref. [14]. We also note that Xm,q diverges as function of the light quark mass at 



t = 0, i.e. at the chiral phase transition temperature. In contrast to the chiral susceptibility, Xm ~ 
dM/dmi, which stays finite in the chiral limit only for t >0, the mixed susceptibility is finite for all 
t ^0. For small values of the light quark mass numerical results for the mixed susceptibilities Xm,q 



may be compared to the right hand side of Eq. (3^). Here all parameters that enter /^(z) are known 
and the only undetermined parameter is Kq. Using a subset of the data samples that have been used 
for the scaling analysis of the order parameter ||7|, ||], we calculated the mixed susceptibility Xm,q 
on lattices with temporal extent = 4 for several values of the quark mass. For this analysis we 
used data sets separated by 50 trajectories. For the lightest quark mass ratio, mi/rris = 1/80, we 
selected 4 and for the three heavier quark mass ratios, mi/nis = 1/10, 1/20, 1/40, we choose 6 
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Figure 2: The mixed light quark number susceptibiUty as a function of the reduced temperature, {T — Tc)/Tc 
(left) and the scaled mixed susceptibility as function of the scaling variable z = t /h^l^^ (right). Shown are 
results obtained at two values of the cut-off, ~A (open symbols) and A^t = 8 (filled symbols), and for 
several values of the light to strange quark mass ratio. On the right hand side, the data is compared to the 
0(2) scaling curve. 



values of the gauge coupling in a narrow temperature interval close to the chiral phase transition 
temperature Tc, i.e. —0.02 < {T — Tc)/Tc < 0.06. Typically this involved about 500 to 950 gauge 
field configurations per parameter set, except for the lightest quark mass ratio where we analyzed 
about 350 gauge field configuration. On each gauge field configuration we calculated the various 
operators necessary to construct Xm,q- 

The calculation of the various operators required inversions of the staggered fermion matrix 
with a large set of random noise vectors. We used 500 noise vectors on each gauge field configura- 
tion and constructed unbiased estimators for the various traces that need to be calculated. All these 
calculations could be performed very efficiently on a GPU cluster 

Results obtained for the mixed light quark number susceptibility, Xm,q^ on lattices with tem- 
poral extent A^^ = 4 are shown in Fig. ^ (left). We clearly see that Xm,q increases in the transition 
region with decreasing values of mi/nis. Using the scaling relation given in Eq. (3.4) we can re- 
scale the data and obtain a unique scaling curve. This scaling curve can be mapped onto the 0(2) 
scaling function /^(z) with a simple multiplicative rescaling factor, 2Kq. The resulting scaling plot 
is shown in Fig. || (right). To check for possible contributions from scaling violating terms we have 
analyzed the data separately for quark mass ratios m//m.s = 1/10, 1 /20 and m//m.s = 1 /40, 1 /80. 
These fits agree within statistical errors. We then determine the curvatures Kq from fits to the 
complete data set. Results of these fits are summarized in Table ^ 

The scaling analysis performed for the mixed susceptibility on lattices with temporal extent 
Nx = A suggests that the determination of the curvature parameter can be reliably performed with 
quark masses m//m4< 1/10. This is in accordance with the scaling analysis of the order parameter 
itself §]. It thus seems to be safe to extract the curvature parameter also at smaller values of 
the lattice spacing, i.e. from our = % data set, by using the smallest quark mass ratio available 
there, mi/nis = 1/20. We have performed calculations at five values of the temperature using 
gauge field configurations on 32^ x 8 lattices generated by the HotQCD collaboration [p^]. For 
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mi/m. 






4 


1/10, 1/20 


0.0598(26) 


3.5 




1/40, 1/80 


0.0573(29) 


1.5 


8 


1/20 


0.0559(35) 


0.4 


4, 8 


all 


0.0591(17) 


2.1 



Table 2: Determination of the curvature of the critical surface of the chiral phase transition in (2 + 1) -flavor 
QCD as function of the light quark chemical potential pL^. The table summarizes fits performed separately 
for two lighter and two heavier quark mass sets as well as the combined data set. 

these parameter sets we have analyzed 300 to 600 gauge field configurations, which were separated 
by 100 trajectories. Again we used 500 noise vectors for the calculation of all relevant operators 
on each of the gauge field configurations. The result of this analysis is shown in Fig. § with filled 
symbols. As can be seen they agree well with results obtained on coarser lattices. 

When rescaling data obtained for Xm,q to the 0(2) scaling curve /^(z) we need to take into 
account errors on the scaling parameters and zo (or ho). This leads to a 10% error for the determi- 
nation of the curvature terms. Performing a combined fit to all results obtained for different quark 
mass values and lattice spacings we obtain Kq = 0.059(2) (4) or equivalently Kb = 0.00656(22) (44). 
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